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Motivated by the inability of Galactic chemical evolution models to reproduce some 
of the observed solar neighbourhood distribution of elements (and isotopes) with atomic 
numbers 6<Z<15, we have revisited the relevant stellar and Galactic models as part of an 
ambitious new program aimed at resolving these long-standing discrepancies. Avoiding 
the use of (traditional) parametric models for low- and intermediate-mass stellar evolution, 
■ we have generated a new, physically self-consistent, suite of stellar models and integrated 
the nucleosynthetic outputs into GEtool, our semi-analytical galactic chemical evolution 
software package. The predicted temporal evolution of several light- and intermediate- 
mass elements (and their isotopes) in the solar neighbourhood - from carbon to phosphorus 
- demonstrate the efficacy of the new yields in reconciling theory and observation. 

1. MOTIVATION 
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The role of asymptotic giant branch (AGB) stars in contributing to the chemical enrich- 
ment of the interstellar medium (ISM) has long been appreciated [dEj- Unfortunately, the 
! dearth of self-consistent grids of non-parametric stellar models for low- and intermediate- 
mass stars (LIMS) ensured that that "appreciation" has remained more qualitative than 
quantitative. Nowhere has this been more problematic than in attempts to understand 
the origin and evolution of the light- and intermediate-mass metals - atomic numbers 
between 6 and 15 - where the discrepency between supernovae yields and observed stellar 
abundances is most acute 

We have recently undertaken an ambitious program of coupled stellar and Galactic 
chemical evolution modeling, preliminary results for which are described here. Of order 50 
stellar models were run using the Mount Stromlo Stellar Structure code; a nucleosynthesis 
post-processing with time- dependent diffusive mixing was then applied, in order to derive 
detailed yield information [HIE]. The parameter space covered was extensive, with masses 
ranging from 1-7 M Q , metallicities ranging from zero to super-solar (for both scaled-solar 
and a-enhanced abundance ratios), as well as varying treatments of mass-loss and reaction 
rates. Elements for which an initial mass function (IMF)-weighted abundance with respect 
to solar varied by more than 0.1 dex from that derived using the canonical LIMS yields 
of [[I] include atomic numbers 6<Z<15 (carbon through phosphorus). 
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These new stellar yields have now been implemented within GEtool [00, a semi- 
analytical galaxy evolution code which treats the formation of the Milky Way within a 
dual-phase infall framework - the first, a rapid infall phase leading to the formation of the 
stellar halo; the second, a more extended phase associated with the formation of the disk. 
Such dual-phase formation scenarios have proven success in reproducing the metallicity 
distrubtion of nearby stars [El El- GEtool includes a sophisticated treatment of chemical 
enrichment in the interstellar medium, with the yields of Type la and II supernovae (SNe), 
Wolf-Rayet stars, and AGB stars, all incorporated. This powerful combination of GEtool 
and the new yields has been demonstrated recently by [ EH HH H2 US] • 

2. RESULTS 

In the following subsections, we highlight just a few specific examples for which the 
new LIMS yields have led to interesting new insights into the distribution of elemental 
and isotopic abundances for atomic numbers 6<Z<15. 

2.1. Carbon 

Over the past 4.5 Gyr, the ratio of 12 C-to- 13 C (by mass) in the local ISM has decreased 
~20%. Attempts at recovering this temporal behaviour in 12 C/ 13 C has led to the sug- 
gestion that significant pollution of the local ISM by novae ejecta must have occurred 
over the lifetime of the Milky Way's disk [El- This conclusion was driven in part by 
the failure of standard LIMS yields to predict the observed decrease in the temporal 
evolution of 12 C/ 13 C over the disk's lifetime - see Fig 1 of [ 14J. 2 Conversely, after in- 
corporating the new grid of LIMS yields into GEtool our standard solar neighbourhood 
model predicts an ~15% decrease in 12 C/ 13 C over the past 4.5 Gyr, in agreement with the 
aforementioned empirical data. Over the lifetime of the disk, the same ratio is predicted 
to have decreased by ~40%. In other words, the new yields appear to obviate the need 
for a putative (significant) nova component to Galactic chemical evolution. 3 

2.2. Fluorine 

The nucleosynthesis pathways for fluorine production in AGB stars involve both the 
helium and combined hydrogen-helium burning phases. The primary uncertainty in the 
net production rate can be traced to the adopted reaction rates of 14 C(a,7) 18 and 
19 F(a,p) 22 Ne, and from the inclusion of partial mixing of protons from the envelope in 
the top layers of the helium intershell region [Q31- Models of Galactic chemical evolu- 
tion which include only 19 F production from neutrino spallation of 20 Ne in Type II SNe 
underproduce the observed [F/O] in the Milky Way and LMC by a factor of ~2 [E?l IT3]. 
After including the new LIMS yields within GEtool, supplemented with the 19 F yields for 
Wolf-Rayet stars at super-solar metallicity [E], we ran our standard Milky Way chemical 
evolution model. The new yields resulted - for the first time - in the successful recovery 

2 We note that the "single star" 12 C/ 13 C models of were scaled upwards by 35% in order to ensure 
an a posteriori agreement with the Sun's 12 C/ 13 C. 

3 The present-day Galactic 12 C/ 13 C gradient is an additional constraint on novae pollution in the ISM. 
We note in passing that our Galactic models, with the new yields, result in a gradient in excellent 
agreement with the Galactic distribution of carbon monoxide [ 1151 . The inclusion of novae (significant 
13 C factories) does not impact significantly on the predicted gradient in 12 C/ 13 C, but the resulting 
zeropoint is approximately a factor of two below that of the present-day Galactic ISM. 
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of the trend of [F/O] observed over ~1 dex in metallicity in the Milky Way and LMC [ 

eng. 

2.3. Neon 

The neon isotope 22 Ne is produced with relative ease in the helium-burning shell of 
AGB stars via the capture of two a-particles onto residual 14 N remaining from the earlier 
operation of the CN-cycle; subsequent thermal dredge-up mixes this 22 Ne into the enve- 
lope. This 22 Ne nucleosynthetic pathway appears to operate independent of metallicity, 
but only over a limited mass range (msi3.0±0.5 M©) [E]. The 22 Ne/ 20 Ne ratio for this 
22 Ne-enriched ejecta (from stars with m~3 M Q and Z<Z & ) is ~3, comparable to the ratio 
expected in Wolf-Rayet stellar ejecta: 22 Ne/ 20 Ne^3-10 for m>40 M and Z>Z Q [Ej. 
The new yields, in combination with our standard Galactic model, result in a factor of ~2 
increase in 22 Ne/ 20 Ne (with no discernible impact on 21 Ne/ 20 Ne) over the past ~7 Gyr 
history of the solar neighbourhood. In contrast, the isotopic ratio 22 Ne/ 20 Ne in Galactic 
cosmic rays has been measured to be anomalously high - ~5x that of the solar wind. 
This observation has been attributed to the cosmic rays having been accelerated from 
superbubbles of metallicity ~3Z© with the accompanying Wolf-Rayet (and Type II SNe) 
ejecta accounting for ~20% of the local ISM (by mass) [HZ]. It would be useful to revisit 
this latter conclusion in light of this previously unappreciated source of 22 Ne - AGB stars 
of mass ~3 M . 

2.4. Sodium 

It is a well-known fact that Galactic chemical evolution models which incorporate 
sodium yields from Type II SNe alone tend to underproduce [Na/Fe] by a factor of ~2-3 
over ~3 dex in metallicity [ 3J. It has also been recognised for some time that the Ne-Na 
chain acting in AGB stars can lead to the production of 23 Na via proton capture on 22 Ne 
[E]. 4 Using our new yields, we have constructed the first Milky Way chemical evolution 
model which includes self-consistently the sodium production from both Type II SNe and 
AGB stars. Regardless of the Type II SNe yields adopted, the inclusion of sodium from 
AGB stars results in fairly uniform 0.2-0.4 dex increase in the predicted [Na/Fe] in the 
ISM for — 2<[Fe/H]<+0, in excellent agreement with observations from [^J^3- 5 

2.5. Magnesium 

While the bulk of magnesium in the Galaxy can be traced to Type II SNe, [ 0] have 
shown that sub-solar metallicity AGB stars can be important contributors of 25 Mg and 
26 Mg isotopes to the ISM. Nucleosynthesis of these isotopes is believed to occur via a- 
capture on 22 Ne triggered by helium shell thermal pulsing. More massive AGB stars 

4 Various nucleosynthesis pathways exist, including the hydrogen- and helium-burning shells, and hot 
bottom burning in massive AGB stars. 

5 It should be noted that a "numerical" (coincidental!) pseudo-degeneracy exists between sodium from 
AGB stars and sodium from stars of mass m=40-100 M Q . For example, as the most massive model 
generated by [|3U] is 40 Mq, if one decides to adopt an upper mass limit of, say, 100 M Q for the IMF, 
one is forced to extrapolate the sodium production from the lower mass models. Because this sodium 
production is a steeply increasing function of stellar mass, linearly extrapolating to 100 M Q can actually 
lead to a predicted [Na/Fe] vs [Fe/H] behaviour which mimics (again, coincidentally) that encountered 
when using the new AGB sodium yields. A careful consideration of very massive star sodium yields must 
be undertaken before we can make any further quantitative statements. 
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(m=4-6 M Q ) may actually burn magnesium via hot bottom burning at the base of the 
convective envelope. Models of Galactic chemical evolution which include only the heavy 
magnesium isotopes returned to the ISM via Type II SNe are significantly discrepant with 
observational data [EJ. Incorporating the Mg isotopic contribution from AGB stars into 
GEtool leads to a factor of 2-3 increase in 25 > 26 Mg/ 24 Mg over ~2 dex in metallicity, in 
agreement with the distribution observed in the solar neighbourhood [I10j. 

2.6. Phosphorus 

Our new stellar models produce 31 P efficiently; when coupled with the phosphorus as- 
sociated with Type II SNe [ |2DJ 1211) an ~0.2 dex enhancement in the predicted solar 
neighbourhood [P/Fe] is seen across ~2 dex in metallicity (— 2<[Fe/H]<+0). While em- 
pirical stellar phosphorus abundances are difficult to determine, a firm upper limit of 
[P/S]<+0 has been placed on the intergalactic medium 6 ~2 Gyr after the Big Bang [ I12j . 
Such an empirical limit was already only 0.3-0.5 dex outside the predictions of chemical 
evolution models generated without an AGB phosphorus component ([ 112] : Fig 6); the 
inclusion of our new 31 P AGB yields into the same damped Lyman-a model means the 
current empirical limit is less than a factor of two outside the model predictions. Future, 
more sensitive, observations of such high-z clouds can therefore, in principle, support or 
refute this specific prediction of our new LIMS models. 
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6 Or at least one damped Lyman-a system. 



